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In-Plane fracture of laminated fiber reinforced composites with 
varying fracture resistance: experimental observations and numerical 

crack propagation simulations 
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Abstract 

A series of experimental results on the in-plane fracture of a fiber reinforced laminated composite 
panel is analyzed using the variational multi-scale cohesive method (VMCM). The VMCM results 
demonstrate the influence of specimen geometry and load distribution on the propagation of large 
scale bridging cracks in the fiber reinforced panel. Experimentally observed variation in fracture re¬ 
sistance is substantiated numerically by comparing the experimental and VMCM load-displacement 
responses of geometrically scaled single edge-notch three point bend (SETB) specimens. The results 
elucidate the size dependence of the traction-separation relationship for this class of materials even 
in moderately large specimens, contrary to the conventional understanding of it being a material 
property. The existence of a “free bridging zone” (different from the conventional “full bridging 
zone”) is recognized, and its influence on the evolving fracture resistance is discussed. The numeri¬ 
cal simulations and ensuing bridging zone evolution analysis demonstrates the versatility of VMCM 
in objectively simulating progressive crack propagation, compared against conventional numerical 
schemes like traditional cohesive zone modeling, which require a priori knowledge of the crack path. 

1 Introduction 

Laminated fiber reinforced composites are finding increased use as structural materials in a variety 
of aerospace and other industrial applications. An appealing property of these composites is their 
high specific strength (strength per unit weight). Even though a strong case can be made for using 
composite structural parts in several areas of a load bearing structure, the requirement to demon¬ 
strate structural integrity and damage tolerance (SIDT) necessitates a fundamental understanding 
of the mechanical response, damage tolerance and damage growth of a load bearing composite struc¬ 
ture. While several investigations have addressed damage growth in the form of delamination crack 
growth, relatively little has been done to understand damage propagation when a crack, or damage 
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in the form of a crack that has severed all laminae of a laminate, is present in a composite structure. 
The influence of such wide area damage on the load bearing ability of a homogeneous and isotropic 
material has received considerable attention in the past; however, a similar effort at understand¬ 
ing issues in a non-homogeneous and macroscopically orthotropic structure is still a problem that 
requires resolution. Because of the different length scales associated with the microstructure of a 
composite material and the resulting composite structure, a multitude of failure mechanisms can be 
simultaneously operative, leading to a very complex damage progression in a composite structure. 
A sharp, through the thickness crack can be present in these composites initially, but, as soon as 
local damage (this can be in the form of matrix micro-cracking) accumulates, crack blunting and 
distributed damage occurs across the highly stressed areas around the initial crack tip. As this 
initial crack starts to grow, a damaged zone of material (bridging zone) evolves in the wake of the 
instantaneous crack tip. Thus, unlike in monolithic materials (such as metals), there is no well 
defined “crack” that can be identified. Instead, a diffused zone of damage is seen to advance. This 
distributed damage results in additional resistance to advancing damage growth, largely contributed 
by fiber bridging and pullout in the crack wake . This enhanced fracture resistance is desirable and 


is a major contributor to the increased toughness of these laminated composites (Cooper (1970), 
Aveston et al.j (19711, Aveston and Kelly (1973), Cox (1991)). 


Analytical models based on linear elastic fracture mechanics (LEFM) have been developed and 
implemented within finite element codes to study a variety of fracture problems. LEFM based ap¬ 
proaches have proven to be effective in predicting crack initiation and subsequent growth in cases 


where material nonlinearity is negligible, and process zones are small (Hertzberg (1983), Xie and 


Biggers (2006a), Xie and Biggers (2006b), Xie et al. (2004), Xie et al. (2005)). However, in het 


erogeneous materials, like laminated fiber composites, the process zone size may be larger than 


any characteristic problem length scale, and thereby, the basic tenets of LEFM cease to hold (Xie 


et al. (2006)). Several mechanisms may contribute to this situation. Micro-cracking, fiber bridging, 
coalescence of voids and other microstructural mechanisms can give rise to a process zone that is 
considerably larger than that permitted for the application of LEFM models. Furthermore, the 
material non-linearity that is induced by these mechanisms leads to a relief of the singular fields 
that would otherwise persist in a strict LEFM setting of an elastic material. A new length scale, 
I*, emerges that is related to a characteristic elastic modulus E, fracture toughness T and cohesive 
strength, CTc, defined as, If I* is larger than any characteristic length scale in the prob¬ 


lem, then, cohesive zone models (CZM) become an indispensable tool for analysis (Pietruszczak and 


Mroz 

(1981), 

Ungsuwarungsri and Knauss 

T^, 

Tvergaard and Hutchinson 

(1992) 

Schellekens 

and Borst ( 

1993), 

Xu and Needleman 

(1994 

, Camacho and Ortiz 

(1996 

»• 


In order to implement a CZM in its simplest form, two fracture parameters are required: a 
fracture toughness (or energy) and a cohesive strength. Both parameters can be experimentally 
determined by coupon level tests, and are subsequently used as material properties for prediction of 
crack growth in other structural configurations. In CZM, an existing crack starts to grow when the 
crack tip stresses reach the cohesive strength, and evolves according to the energy available in the 
deformed system to create additional cracked surface. Based on these two parameters, a cohesive 
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traction-separation law is assumed in numerical simulations of crack growth. This micro-mechanical 
law relates the evolution of crack face tractions with the crack face opening displacement, and nu¬ 
merically manifests the resistance offered to crack advancement within the cohesive zone. This two 
parameter dependent evolution, unlike the solely fracture toughness based LEFM schemes, better 
represents the physics of crack growth in materials with significant process zones, and it has been 
widely used in many numerical implementations, like the discrete cohesive zone method (DCZM) 


(Xie and Waas (2006), Xie et al. (20061). However, all implementations of CZM methodology require 


a priori knowledge of the “intended” crack path. This prevents CZM from being applied to a wide 
range of problems involving arbitrary crack path evolution. Conventional finite element based imple¬ 


mentations of CZM (Nguyen et al. (2001)) constrain the crack path evolution to element interfaces, 
as standard finite elements lack the ability to represent cracks within the element domain. This 
restriction places a limitation on those problems for which the crack path direction is not known a 
priori. 


In addition to the above drawbacks of CZM models, recent work by Jin and Sun (2005), Jin and 


Sun (2006), has addressed some basic issues pertaining to the reconciliation of LEFM and CZM. In 


particular, two issues have been considered; (1) when CZM is used, the placement of CZM elements 
along the intended crack path can lead to an alteration of the stiffness of the original body that is 
to be studied, and, (2) The traction-separation laws used for traditional CZM modeling, which start 
with a vanishing traction at vanishing separation, may be in conflict with the presence of an intense 
stress field that was present in the original body that is being modeled. 


To order to address these basic issues and to circumvent the numerical restrictions on unrestricted 
and objective simulation of crack propagation in materials, a micro-mechanics based, mesh indepen¬ 
dent numerical technique for simulating crack propagation is essential. Standard finite elements fail 
to accomplish this task as they lack the ability to capture the discontinuous displacement modes 
involved in crack propagation problems. However, in recent years finite elements with discontinuities 
(enhanced finite elements) have gained increasing interest in modeling material failure, due to their 
ability to capture the specific kinematics of a displacement discontinuity(like cracks) through addi¬ 
tional discontinuous deformation modes. In discontinuous displacement enhanced finite elements, 
the crack path is present inside the elements, unlike cohesive zone methods which are restricted to 
crack propagation along element interfaces. The ability of the enhanced finite element to encompass 
a crack path, leads to objective simulation of crack propagation without mesh bias. Depending on 
the support of the enriching discontinuous displacement modes, the enhanced finite elements are 


popularly classified as element enrichment methods 

Armero and Garikipati (1996), Garikipati and 

Hughes 

1998), 

Jirasek (2000 

, Borja and Regueiro 

(2000), Oliver and Huespe (2004|, Mosler and 

Meschke 

(2004) 

, Gasser and Holzapfel (2003)) and nodal enrichment methods (XFEM, jBelytschko 

et al. (2001), Belytschko et al. 

(2003), Moes et al. (2000), Wells and Sluys (2001)). Interested readers 

are referred to ( 

Jliver et al. (2006), for detailed discussion and comparison of these methods. Though 


these enhanced methods provide a general numerical framework for simulation of crack evolution, 
the actual micromechanics implementation which incorporates the physics of crack formation is 
wide open. In this context, we present the Variation Multiscale Cohesive Method (VMCM), which 
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is an enhanced finite element method containing elemental displacement field enrichment, naturally 
arising out of the variational multiscale formulation presented in Garikipati and Hughes (1998), 


Garikipati (2002), seamlessly embedding the cohesive nature of crack path evolution. In this paper. 


the VMGM method advanced by the authors is briefly presented, and it is used to study through 
the thickness crack propagation in fiber reinforced laminated panels. A more detailed presentation 
of VMGM is available in other related studies by the authors (Garikipati (2002[), [Rudraraju et al. 


(2008)), where the method is referred to as the variational multiscale method (VMM). 


The paper is organized as follows: In Section the application of VMGM to study crack propa¬ 
gation is motivated through a discussion of experimental results in scaled SETB specimens, loaded 
to failure. Section [^presents the details of VMGM implementation using the finite element method. 
Simulations of SETB tests with varying specimen size are provided in Section Comparison of 
experimental results, against the predictions of VMGM are discussed in Section while concluding 
remarks are presented in Section 


2 Motivation 

The primary motivation for the current work is the experimentally observed scaling in the in-plane 
fracture resistance^ with increase in specimen size of single edge notch three point bend (SETB) 
specimen, subjected to primarily Mode-I crack tip conditions. The material used in all the experi¬ 
ments herein is a carbon fiber/epoxy [—45/0/ -f 45/90]6s laminated fiber reinforced composite with 
a volume fraction of 0.55, and whose lamina and laminate properties are given in Table 

The SETB configuration used in this study is shown in Figure]^ SETB specimens were cut from 
the composite panels using water jet facility. The notch was introduced and a knife edge was used to 
introduce a sharp starter crack. The specimens were supported on rubber rollers both at the loading 
and support points to minimize any local inelastic deformation. The specimens were loaded on a 
specially designed loading frame with anti buckling guide rods that prevents out of plane movement 
of the specimens. The specimens were loaded at a rate of O.Olmm/sec using hydraulically operated 
MTS testing machine and were loaded until failure. Load was measured by a load cell and the load 
point displacement was measured in between the top and bottom loading rollers using an LVDT. 

Five specimen sizes with geometrically scaled planar geometry and fixed thickness were considered. 
Multiple specimens of each size were tested to significantly capture the failure response envelope. 

The load - load point displacement (PA) responses of these specimens are shown in Figure]^ The 
apparent in-plane fracture resistance (hereafter referred to as the “in-plane fracture resistance”, or 
simply as “fracture resistance”) of each specimen is determined by inverse modeling, i.e, by de- 

Traditionally, fracture resistance is defined in terms of the energy released per unit area of crack surface formation. 
But, unlike in most monolithic materials, there is no well-defined “crack” that can be identified in a laminated fiber 
composite panel, where a crack like feature propagates severing all laminae of the laminate. Then, by unit area of crack 
surface, we mean unit area of completely failed (both matrix and fiber failure) surface along the crack like diffused zone of 
damage. The usage of the term “fracture toughness” is resisted, to avoid confusion with the traditional LEFM approach, 
where it is often regarded as a material property. 
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termining the value of fracture resistance input required in the VMCM model (Section 4.1), such 
that the experimental and numerical PA curves and the corresponding instantaneous crack lengths 
nearly match. The multiple tests for each SETB size show some scatter in the PA response. Hence, 
an averaged fracture resistance, 3lavg, of each specimen is calculated by averaging the individual 
in-plane fracture resistance values obtained by inverse-modeling from each of the PA curves of that 
specimen ' * As will be demonstrated and discussed in Section ^avg captures the correct failure 
response, and it is also much easier to determine. 


Table summarizes the observed scaling in the PA response, and hence in the value of. 


These quantities P*, A* and 


-aug 


are fixed reference values, of which , 


•'avg 


r-avg- 

is the value of fracture 


resistance obtained in similar tests on standard Compact Tension (CT) specimens (Figure 3(a)). 
The scaling observed in ^avg is very significant, because in all the specimen sizes considered, it has 
been experimentally observed that there is formation of a Full Bridging Zone (FBZ) (Figure]^, and 
the length of FBZ scales up with specimen geometry. FBZ is assumed to have formed when fibers 
at the location of the initial crack tip have failed/completely pulled out, leading to zero tractions 
across the crack faces at the initial crack tip (Cox ( 1991[ )). FBZ formation is also referred to as the 
formation of a “stable process zone”. When the fibers at the initial crack tip have not yet failed, 
the bridging zone is still evolving, and herein we refer to this as the partial bridging zone. 


The formation of a FBZ in all specimen sizes considered and the observed scaling in its length 
are contrary to the usual expectation, wherein, there is no FBZ formation prior to complete failure 
in smaller specimens, and when FBZ formation does occur in large enough specimens, the FBZ 
length is a material property. This assumption of a constant length FBZ formation can lead to 
incorrect predictions of strength and reliability, because in some classes of materials, like the one 
being considered here, a FBZ formation can be misinterpreted as the formation of a “converged 
process zone”, whereas, in reality, larger size specimens might still lead to longer process zones, and 
hence greater FBZ length. Here, the phrase “converged process zone” refers to a process zone that 
does not increase in length if the specimen size is increased beyond a certain size. That is, the 
converged process zone length is invariant with respect to further increase in specimen dimensions. 
In order to substantiate and understand the mechanics behind such scaling behavior and FBZ for¬ 
mation, VMCM based simulations were conducted. Before proceeding to the numerical simulations, 
the VMCM framework is briefly presented in the following section. 


^It is interesting to note that in this material, the value of fracture resistance obtained by inverse modeling nearly 
matches the value obtained by normalizing the total area under the PA curve by the area (thickness times total crack 
length) of the crack surface formed. 

^One may suggest the R-curve approach of LEFM as a better alternative. However the failure response of this 
material is stochastic, and lacks post-peak PA response repeatability, which is important for obtaining consistent R- 
curves. Averaging R-curves across specimens of each size is possible, but it is of little advantage over ^avg- Furthermore, 
the total energy released may include dissipation due to other mechanisms, like crack tip inelasticity, but for this class of 
materials, the crack formation energy, which includes fiber pullout, fiber breaking and matrix cracking, is the dominant 
component of the total energy dissipated. 
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3 Mathematical Formulation 


The standard weak form of the balance of linear momentum over 


the domain (Figure]^ is given 


by, 


: (T dV = / w- b dV + / 


w- T dS. 


( 1 ) 


JQ JQ J 

where cr is the stress, w is an admissible displacement variation, is the symmetric gradient of 
the variation, T is the external traction and b is the body force. 


In the standard finite element formulation of continuum mechanics, the displacements are at 
least (7° continuous. But in a wide class of problems (shear banding, fiber kink banding, transverse 
crack formation, delamination initiation are some examples), the displacement field can be discon¬ 
tinuous. In such cases, the displacement field can be decomposed into continuous coarse scale and 
discontinuous fine scale components (Figure]^. Such a decomposition is also imposed upon the 
displacement variation, w. The decomposition is made precise by requiring that the fine scales, u' 
and w', vanish outside of some region Id', which will be referred to as the microstructural or fine 
scale subdomain. This decomposition is written as. 


U = "b 

coarse scale fine scale 

The corresponding scale separation in the displacement variation is given by. 


( 2 ) 


W = ”b 

coarse scale fine scale 
u', w'eS'— {v'\v'— 0 on ^\int{n')} 

Substituting the above decomposition into 0. and using standard arguments, the weak form 
can be split into two separate weak forms. One, involving the coarse scale variation, w, and the 
other, involving only the fine scale variation, w'. 


; (T dV = / w- b dV -I- 


w- T dS. 


( 4 ) 


Id fit 


/ V"u)':crdV=/ w' bdV+ w'-T dS. (5) 

Jfl' Jfl' Jdfl't 

This procedure results in the fine scale weak form ([^, defined only over fl' (Figure]^. This result 
is crucial since it lends itself naturally to the application of desired micromechanical descriptions 
restricted to the microstructural region, O', and not the entire domain O. The scale separation in 
u is contained in cr = C: (V^ii -I- V®m'), where C is the elastic stiffness tensor. 


We wish to use an appropriate micromechanical law by which the fine scale solution, tt', can be 
expressed in terms of u and other fields in the problem. Below, we will show how such a microme¬ 
chanical law can be embedded into the formulation using the weak form 0. The final step involves 
elimination of the fine scale displacement, it', from the problem by substituting its relation to u in 
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the coarse scale weak form Q. Thus, the fine scale solution does not appear explicitly; however, its 
effect is fully embedded in the resultant modified weak form. 


We choose to contain the crack surface T on which u' is discontinuous. Invoking standard 
variational arguments, the weak form of the fine scale problem can be reduced to the following 
statement of traction continuity (Garikipati (2002)): 


|crn]p = 0 (6) 

where |.] is the discontinuity in the quantity and n is the normal to the crack surface, T. Writing 
the traction on T in terms of components T„ and Tm along n and m respectively (Figure]^, the 
traction continuity condition can be expressed as. 


T„n + T^m = (Tn|r- (7) 

The traction cm|p-, is determined by the macromechanical continuum formulation. The evolu¬ 
tion of Tji and Tm is governed by the micromechanical surface law (Section |^A)). This law, which 
emerges at a finite value of traction, is related to the displacement discontinuity which is the sep¬ 
aration between the surfaces. There are two traction laws across a planar (2D, line) surface and 
the displacement discontinuity |m] can be expressed in terms of the normal opening, |M].n, and 
tangential slip, litj.m, across T. 


We now consider a specific functional form for the micromechanical model, emerging at a non¬ 
vanishing traction. 


Tn = Tno-'Hnluj.n, - "HmM.m (8) 

where and Tm^ are the maximum values of T„ and Tm admissible on T (Figure]^, |M].n > 0 is 
the normal jump (Mode-I type crack opening) and litj.m is the tangential slip (Mode-II type crack 
face slip) along the elemental crack face, Hn and Tim are the softening moduli for the Mode-I and 
Mode-II crack opening evolution, respectively. Consistency between the micromechanical law and 
the macromechanical continuum description is enforced by ([^ via ([^. 

Substituting (|^ in ([^ and dispensing with the explicit indication of crn|p-, 


{Tno - 'Hn{lul.n))n + {Tmo - 'Hm{lu\-m))m - crn = 0 
Expanding ([^ up to first order terms, in order to solve for u': 


(9) 


{Tna - 'Hnilu\-n))n T {Tmo - 'Hm{lul.Tn))Tn - (TTl 
-'H„(5|M].n)n - nm{Slul.m)m 

-{C: {VSu + VSu'))n = 0 (10) 
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where the first line in (101 represents a zeroth-order approximation to (|^ , and the remaining terms 
are the first order corrections. Using u' = |m] Cp, where Cp is the fine scale interpolation (Figurej^, 


converts (10) into a linear equation in ^|it] which can be solved, and then the incremental fine scale 
field is obtained from du' = <5|tt]Cr. Formally, it is represented as. 


— .^[1/., (T, ^n: ^m] 


( 11 ) 


Extending the incremental formulation to cr, which in a general nonlinear problem can be ex¬ 
panded up to first order as cr = cr^^^ + C: (V^m -b where cr(°) is the converged value of 

cr in the previous solution increment, and substituting u' in Q, we obtain the coarse held weak 
form which is independent of the hne scale displacement u'. On solving for 6u, the incremental 


hne scale held 6u' can be recovered via (11). Iterations are to be performed: + 5u, 

= -b 5u', until a converged solution is obtained. From (^4L(10) and (11), it should 


u 


be clear that the VMCM method results in an embedding of the micromechanical surface law into 


the coarse scale weak formulation. Interested readers are referred to Garikipati (2002) for a more 
detailed discussion of the numerical framework. 


Having briehy presented the numerical formulation, we direct our attention to simulations, and 
the mechanisms involved in the in-plane fracture of laminated hber composites. 

4 Numerical Simulations 


4.1 Micromechanical surface law 

The micromechanics of crack propagation is embedded into the macroscopic formulation of the 
VMCM by enforcing the linear traction evolution relations (©• Due to the stochastic behavior 
observed in experiments and nonlinearity of response (Section [^ , it is difficult to obtain precise 
functional forms for "Hn and T-Lm- The accepted methodology in the cohesive zone community 


(Hillerborg et al. (1976), Dugdale (1960), Barenblatt (1962)) is to assume a suitable functional form 
of the traction evolution response, referred to as the traction separation law, anchored by experi¬ 
mentally determined values of r„g, Mode-I I^avg and Mode-II ^avg- In the current study, a 
linear traction separation law (Figure]^ is assumed for the evolution of There is no theoretical 
or physical rationale in assuming a linear traction separation law as against any other functional 
forn0 The actual shape of the traction separation law is still a topic of active research (Scheider 
[and Brocks (2003)) and some methods have been devised to measure the law for a limited class of 


configurations (Andersson and Stigh (2004), Sorensen and Jacobsen (1998)). 


order to avoid additional physical complexities arising from mode mixity conditions only test cases undergoing 
pure Mode-I evolution are considered in the current study, and hence there is no evolution of Tm- However, the VMCM 
framework is general and works just as well for problems involving curved crack propagation and mixed-mode evolution 


(Garikipati (2002l). 

^For this class oi 
tractions separation curve. 


^For this class of materials the simulations are more sensitive to the values of Tng and Siavg, than the shape of the 
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The value of T„j, is obtained by experiments on double notched tension specimens (Figure pl(b)[ ). 
This configuration was selected because the stress state across the entire crack face is almost uniform 
and specimen failure is instantaneous. Thus, the critical load divided by the total crack area gives 
a fairly accurate estimate of the critical traction across the crack faces. It was observed that this 
value was independent of the specimen size. The ^Havg values are given in Table ([^. With these two 
inputs, the five SETB specimen sizes given in Figure are simulated using the VMCM. Further, 
the crack is propagated perpendicular to the maximum principal stress direction. 


4.2 Simulations 

The criterion selected in this work for validating the VMCM methodology is the comparison of 
VMCM and experimental PA responses of SETB specimens (Figure]^. This criterion is chosen 
because the PA response reflects the macroscopic response of the structure to external loads, and 
often in structural design the peak load is the value of primary interest. A representative finite ele¬ 
ment simulation mesh is shown in Figurej^ Figurej^shows the PA curves extracted from VMCM 
simulations of SETB specimens with a constant, 3iavg = input for all sizes considered. As 

seen in the figure, the load displacement response is significantly captured for smaller specimens 
(because is close to ^avg value for Size-1 and Size-2 specimens), but the peak load is severely 
under-predicted for larger specimens. This is expected, as the experimental values of ^^avg for large 
specimens are much higher then the input value (Tablej^. Now the simulations are conducted 
with the experimentally measured ^avg values given in Table and the results are plotted in Fig¬ 
ure [T^l 

As seen from Figure [T^ the VMCM simulations accurately reproduce the macroscopic response 
of the SETB specimens when appropriate ^avg values are used as inputThis demonstrates that: 

(1) The VMCM methodology has the ability to numerically simulate progressive damage propaga¬ 
tion, and the mechanics of bridged crack evolution. (2) In spite of a multitude of failure mechanisms 
operating simultaneously, leading to a very complex evolution of fracture resistance, the single val¬ 
ued estimate, ^avg, of the fracture resistance is appropriate for numerical simulations, at least in 
this class of materials!!] 

Having presented the relevance of VMCM and ^avg, we now turn to the actual mechanics of 
the observed scaling behavior: Why does the fracture resistance (and hence ^avg) increase with 
specimen size?. As a start, we use VMCM simulations to study the effect of high fracture resistance 
on the PA response of this material model. Figure [T3| shows the effect of varying ^avg on the peak 
PA response of the specimens. As seen, increase in ^%avg leads to plateauing of the peak response. 

The peak load increases only slightly, but the tendency for softening past the peak is suppressed. 

®The small pre-peak stiffness difference in the experimental curves of some specimens is due to some amount of material 
crushing occurring at the loading points, and the post-peak variations are an artifact of the stochastic nature of failure 
in these materials. 

^Even though the linear traction separation law simplifies a set of complex failure mechanisms, the dissipated energy 
is well-represented. This quantity appears to control the fit with experimental PA curves. 


9 





One inference from this observation is that in this class of materials, high fracture resistance has 
only limited influence in increasing the peak load carrying ability of the specimens. Beyond a certain 
range of fracture resistance, the specimen geometry starts to play a significant role. This is evident 
in Figure |14[ where the bridging zone formation and movement is depicted. As seen in the figure, 
after the formation of a FBZ, the zone propagates for some distance with almost a constant length. 
But as the FBZ approaches the end of the crack path, the compressive stresses encountered by the 
growing crack, caused due to bending of the SETB specimen about the loading point, inhibit further 
FBZ propagation, and hence the FBZ length decreases]^ 


Based on the simulation results depicting the influence of 3iavg on the PA response (Fig¬ 
ure 11 13) and the FBZ formation (Figure [I^, we present the following assessment of the failure 
behavior of this class of materials. It appears that this class of material is capable of achieving a very 
large bridging zone length, when there is no geometric constraint (both, ahead of the crack tip, and 
in the crack wake undergoing fiber pullout) on its natural bridging zone formation tendency. The 
corresponding maximum “full bridging zone” length is referred to as the “free bridging zone” (FrBZ) 
length, which may be considered a material property. In this FrBZ, a multitude of mechanisms (fiber 
failure, debounding, pullout) are active over large lengths, leading to a high fracture resistance. But 
in specimens with sizes comparable to these large lengths, the FrBZ formation mechanisms are in¬ 
hibited and the specimen fails to manifest the FrBZ, but instead forms a smaller FBZ. The smaller 
the specimen size, larger the inhibition to free bridging zone formation, and thereby, smaller the 
corresponding FBZ length manifested. This reduced FBZ length leads to lower fracture resistance, 
and hence lower peak load. 


5 Discussion 


The material in consideration is composed of a quasi-brittle matrix and reinforcing carbon fibers. 
For this material, experiments show at first a matrix crack forming, leading to enhanced loading 
on the adjoining fibers, which suffer greater strain until they reach their failure stress. Ultimate 
failure occurs at weak points either at the matrix crack interface or inside the surrounding matrix. 
Fiber failure at the matrix crack interface leads to complete loss of load bearing capacity at the 
respective fiber. If the fiber fails at a weak point inside the intact matrix, it leads to hber pull out, 
which involves debonding failure and frictional resistance to the fiber movement (Cooper (1970)). 
As stated earlier, the process zone over which these mechanisms occur is termed the FBZ. These 
additional dissipation mechanisms contribute to an increase in fracture resistance until the fiber 
is completely pulled out of the matrix. These mechanisms are schematically depicted in Figure 
Theoretical models have been developed to study fiber failure and pullout, and its effect on fracture 
resistance (Cooper (1970), Aveston et al. (1971), Aveston and Kelly (1973), Cox (1991)). These 


®The maximum FBZ length corresponds to the peak load. Any further progress of the FBZ has to be accompanied 
by a reduction in the FBZ length, and hence drop in P, to satisfy force and moment equilibrium of the specimen. Any 
further reduction in FBZ length due to compressive stresses ahead of the crack tip, is in addition to the reduction of the 
FBZ length necessitated due to equilibrium conditions. 
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models are useful in gaining insight to the individual contributions of various constituent mechanisms 
towards fracture resistance, but in a real material the observed fracture resistance is a combination 
of these processes. Furthermore, the actual crack path formation is significantly stochastic, and 
it is prohibitively difficult, if not impossible, to delineate the contribution of each process. Thus, 
a macroscopic representation of fracture resistance, as done here, provides a convenient means to 
model through the thickness crack growth in this class of materials. 


In Section 4.2 the effect of varying ^%avg on the peak load of the specimen was briefly discussed. 
Figure 13 depicts the limited influence of 3^avg in increasing the peak load of the specimen. As 
evident from the figure, beyond a certain range, increase in 3iavg causes a plateauing of the PA 
response for the following reasons: 

(1) Another parameter, (Equation]^ Figure]^, also influences the peak load magnitude. For 
any given FBZ length, the load carrying ability of the FBZ may be determined by adding the in¬ 
stantaneous load carried by all the individual fibers making up the FBZ. But the instantaneous load 
(Tn) carried by each fiber is dictated by Equation]^ and thus influences the peak load of the 
specimen. 


(2) Though 3iavg and determine the maximum load carrying ability of the material, the actual 
evolution of the PA response has to be in accordance with the force and moment equilibrium con¬ 
ditions. 


It is instructive to place the results that have been presented in the context of energetic size effect 
laws that have been developed by Bazant and Kazemi (1990). When geometrically scaled specimens 
are stressed to failure, it has been shown that many types of quasi-brittle materials obey a Type 
2 size effect law that is described by Bazant and Kazemi (1990). As shown there, performing a 
regression analysis of just the maximum loads for the SETB specimens (Figure]^ leads to estimates 
of FrBZ length and S^max in accordance with the scaling law. Using this ^%raax value as input, 
the VMCM simulations predict the FrBZ length which is in good agreement with the scaling law 
estimate (FrBZvMCM/FrBZ5coHn3Laiu=l-08). This comparison adds further confidence and validity 
to the findings presented here through the use of the VMCM method. 


6 Conclusions 

In this paper, the VMCM method has been used to study in-plane fracture and damage growth 
in laminated fiber reinforced composite panels. Motivated by a series of experimental results, a 
novel approach that circumvents the drawbacks of traditional cohesive zone modeling approaches 
has been developed. The method has the advantage of being able to predict non-self similar crack 
growth along paths that do not need to coincide with element edges. Furthermore, knowledge of the 
intended crack growth path is not needed, and in fact is an outcome of the VMCM method. The 
results presented and the ensuing discussion demonstrate that: 

(1) In-plane fracture of fiber reinforced laminates, unlike in metals, exhibits geometry and loading 
dependent fracture resistance. 

(2) Averaged fracture resistance {^avg) values reflect the macroscopic fracture resistance of these 
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specimens. 

(3) The existence of a “Free Bridging Zone”(FrBZ) and maximum fracture resistance {^max) are 
characteristics of this material. The geometric inhibition of FrBZ evolution and associated conver¬ 
gence of ^avg to ^max lead to the observed scaling behavior. Not accounting for the existence of 
FrBZ and ^max posses the risk of incorrect predictions of strength and reliability, because in some 
classes of materials, like the one being considered here, the FBZ formation may wrongly signal the 
formation of a converged process zone (FrBZ). 
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Table 1: Lamina and laminate properties of carbon fiber/epoxy [—45/0/ + 45/90]gs laminated fiber 
reinforced composite. 


Laminate 

Lamina 

Exx- 51.5 GPa 
Eyyi 51.5 GPa 
Gxy: 19.4 GPa 
Uxy: 0.32 

E\\'. 141 GPa 

E 22 -. 6.7 GPa 
Gi 2 -. 3.2 GPa 

1^12: 0.33 


Table 2: Scaling observed in the SETB specimen experiments. 


Size 

Geometry scaling 
(Figure p]) 

Peak load 
P/P* 

Load point displacement 
A/A* 

Fracture resistance 

^avg/^avg 

1 

1 

0.27 

0.1 

1.08 

2 

1.5 

0.4 

0.15 

1.23 

3 

2 

0.6 

0.2 

1.84 

4 

3 

0.81 

0.28 

2.46 

5 

4 

1.0 

0.37 

2.58 


|P,A 



Size 1: X (as shown) 
Size 2: 1.5X 
Size 3: 2X 
Size 4: 3X 
Sizes: 4X 


Figure 1: Single Edge Notch Bending(SETB) specimen configuration used for validating VMCM sim¬ 
ulation results. Size 1 has the dimensions shown in figure, other sizes are scaled versions of this base 
size. All specimens have a nominal thickness of 6.35mm. 
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Crack opening in Size 5 specimen Crack opening in Size 1 specimen 



Fiber bridging in Size 5 specimen 


Figure 2: Comparison of crack tip opening in failed SETB size-1 and size-5 specimens. 



P,A 


in 

ai 


P 


6.3 mm 


1 mrn 
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(b) 


Figure 3: (a) Compact tension specimen (CTS) configuration used to obtain ^avg for Mode I crack 
propagation, (b) Double notch tension specimen configuration used to obtain critical cohesive stress 
value. All specimens have a nominal thickness of 6.35mm. 
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Figure 4: Experimental Load (P) - Load point displacement (A) curves obtained for various sizes of 
SETB specimens. Multiple specimens of each size were tested to capture the envelope of the failure 
response. P* and A* are fixed reference values. 



Eigure 5; Schematic of scale separation, u is the coarse scale displacement field and u' is the local fine 
scale enhancement. 



Eigure 6: Decomposition of continuum body into region where coarse scale and fine scale displacements 
are defined. D is the domain of the problem, F is the displacement discontinuity (crack), D' is the 
support for the displacement discontinuity and n, m are the normal and tangent to F respectively. 
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h h h 

Figure 7: Discontinuous shape function used to resolve the displacement jump shown in a one¬ 
dimensional setting. It is constructed by superimposing a discontinuous function, Hr, on a regular 
polynomial, N. h is the element dimension. 



Figure 8: Linear micro-mechanical surface law for normal crack opening 





Figure 9: Stages involved in fiber composite cracking. A: Matrix cracking - Fiber bridging, B: FBZ 
formation, C: FBZ propagation. 
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(a) Original mesh 





(b) Deformed mesh: The red line represents the crack path 


Figure 10: SETB specimen mesh, with crack tip elements magnified in inset. 
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Figure 11: Load-Displacement {PA) response, extracted from VMCM simulations with a constant, ^avg — ^avg^ input for Size 
1-5 SETB specimens, compared with the experimental curves. The P and A values have been normalized with fixed reference 
values. 
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Figure 13: Effect of increasing fracture resistance input on the Load-Displacement (FA) response of VMCM simulations, ^avg i® 
a constant reference value. Also plotted is the corresponding FA response obtained with the experimentally determined fracture 
resistance value. The F and A values have been normalized with fixed reference values. 
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stabilization and movement. 




























































